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SUMMARY  OF  PROGRESS 


The  long  term  goal  of  this  research  is  to  develop  a  thorough  understanding  of 
thermoacoustic  phenomena.  During  FY  1991,  our  investigations  concentrated  on 
five  areas:  an  experimental  study  of  edge  effects  in  thermoacoustic  couple 
measurements,  the  application  of  porous  media  techniques  to  thermoacoustic 
prime  movers,  development  of  a  non-boundary  layer,  non-short  stack  model  for  the 
power  output  of  a  thermoacouistic  prime  mover,  the  study  of  the  transition  to  steady 
state  oscillation  in  a  thermoacoustic  prime  mover  above  onset  of  self-oscillation, 
and  the  study  of  energy  distribution  and  dissipation  in  finite  amplitude  standing 
waves.  The  first  two  projects  straddled  FYs  1990  and  1991.  The  remaining  three 
were  initiated  in  FY  1991.  A  brief  discussion  of  these  investigations  is  given  in  this 
report. 


1.  Experimental  Study  of  Edge  Effects  in  Thermoacoustic  Couple  Measurements 

In  previous  years  we  devoted  a  large  effort  to  measuring  the  temperature 
difference  developed  across  a  simple  class  of  thermoacoustic  engine,  called  a 
thermoacoustic  couple  or  TAC.  The  results  of  these  extensive  measurements  are 
discussed  in  detail  in  Ref.  1.  To  summarize  the  findings,  large  discrepancies 
between  the  predicted  and  measured  temperature  differences  app>e8U’  jI  high  acoustic 
amplitudes.  Although  the  reasonfs)  for  these  discrepancies  is  (are)  unknown,  some 
possible  explanations  are  suggested  in  Ref.  1.  One  possibility  involves  the  TAC 
design.  The  temperature  difference  was  measured  with  a  series  of  thermocouples 
(called  a  thermopile)  placed  along  the  front  and  back  edges  of  the  TAC.  Because  of 
its  proximity  to  the  edge,  the  thermopile  may  be  sensitive  to  effects  that,  though 
perhaps  causing  local  deviations  in  the  temperature  profile,  do  not  affect  the 


temperature  profile  in  interior  regions  of  the  plate.  In  other  words,  although 
considerable  deviations  are  observed  between  theory  and  measurements  made  at 
the  plate  edges,  the  interior  (and  majority)  of  the  plate  may  behave  according  to 
predictions. 

To  investigate  whether  edge  effects  are  the  cause  of  the  discrepancies,  we  have 
constructed  a  TAC  with  two  thermopiles,  whose  junctions  do  not  lie  along  the  edge, 
and  made  the  same  type  of  measurements  as  described  in  Ref.  1.  These 
measurements  are  discussed  in  detail  in  Refs.  2  and  3.  The  results  show  that  the 
temperature  proHle  in  the  interior  of  the  TAC  behaves  the  same  as  that  measured  at 
the  edge.  Therefore,  the  irregularities  previously  observed  are  not  an  artifact  of 
placing  the  thermopile  along  the  edge.  The  irregularities  may  still  be  the  result  of 
the  edge  (e.g.,  turbulence  generated  at  the  edge),  but  they  are  not  isolated  at  the  edge. 
They  extend  at  least  half  way  to  the  center  of  the  TAC. 

Although  the  cause  of  the  discrepancy  between  theory  and  measurement  is  still 
unknown,  the  number  of  possibilities  has  at  least  been  reduced.  Of  the  remaining 
possible  causes,  turbulence  and  acoustic  streaming  are  the  most  probable  candidates. 


2.  Application  of  Porous  Media  Techniques  to  Thermoacoustic  Prime  Movers 

In  FY  1990,  we  began  to  investigate  thermoacoustic  prime  movers  below  the 
onset  of  self-oscillation.  The  quality  factor  Q  of  a  helium  filled  prime  mover  was 
determined  from  measurements  of  its  frequency  response.  Measurements  were 
made  for  mean  pressures  ranging  from  170  -  500  kPa  and  temperature  differences 
across  the  stack  ranging  from  zero  to  that  required  for  onset  of  self-osdllation.  As 
the  temperature  difference  is  increased,  the  Q  increases  indicating  a  decrease  in 
attenuation  across  the  stack.  Analysis  of  the  data  and  refinements  to  the  theory 
continued  in  FY  1991  as  discussed  below,  which  is  a  summary  of  a  complete 


disoission  presented  in  Ref.  4. 

The  results  of  the  measurements  are  explained  in  terms  of  a 
counterpropagating,  plane  wave  analysis,  based  on  techniques  commonly  used  in 
porous  media  investigations.  The  results  of  the  calculations  are  shown  in  Figs  1-3, 
along  with  the  experimental  data.  These  figures  show  graphs  of  1/Q,  which  is 
proportional  to  the  net  attenuation  in  the  prime  mover,  versus  temperature 
difference.  The  figures  correspond  to  mean  gas  pressures  of  170,  376,  and  500  kPa, 
respectively.  The  individual  symbols  represent  the  measurements  and  the  lines 
show  the  results  of  the  calculations.  The  agreement  between  the  theory  and 
experiment  at  the  two  higher  pressures  is  very  good.  There  are  some  noticeable 
differences  at  the  lowest  pressure.  The  theory  over  predicts  the  attenuation  at  low 
temperature  differences  and  imder  predicts  the  onset  temperature  (1/Q  =  0). 

Measurements  were  also  made  for  the  second  and  third  longitudinal  modes  of 
the  prime  mover  at  a  mean  gas  pressure  of  170  kPa.  The  data  and  calculations  are 
shown  in  Fig.  4.  The  x’s  and  solid  line  correspond  to  the  second  mode,  while  the 
open  squares  and  dashed  line  correspond  to  the  third.  It  was  not  possible  to  make 
measurements  near  onset  for  these  modes,  because  the  onset  temperatures  exceed 
that  of  the  fundamental.  The  overall  agreement  is  good.  The  tendency  to  over 
predict  the  attenuation  at  zero  temperature  difference  decreases.  No  conclusions 
can  be  drawn  concerning  the  ability  to  predict  the  onset  temperature. 

The  theory  presented  in  Ref.  4  b  intended  only  to  be  approximate  in  that  it  does 
not  consider  the  complications  of  the  temperature  dependence  of  thermophysical 
properties  of  the  prime  mover  components  in  a  detailed  manner.  However,  it  does 
show  that  a  counterpropagating  wave  analysis,  commonly  used  in  porous  media 
research,  does  offer  a  relatively  good  explanation  of  the  data.  Collaborators  at  the 
University  of  Mississippi,  have  developed  a  more  thorough  theoretical 
development^  and  have  used  our  data  to  verify  it.^ 
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Figure  1.  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover 
stack.  The  symbols  represent  the  data,  while  the  line  represents  the  results  of  the 
calculations.  The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  170  kPa. 


figure  2.  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover 
stack.  The  symbols  represent  the  data,  while  the  line  represents  the  results  of  the 
calculations.  The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  376  kPa. 
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Figure  3.  Graph  of  1  /Q  vs  the  temperature  difference  across  the  prime  mover 
stack.  The  symbols  represent  the  data,  while  the  line  represents  the  results  of  the 
calculations.  The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  500  kPa. 


Figure  4.  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover 

stack  for  the  second  and  third  longitudinal  modes.  The  x's  and  solid  line 
cmrespond  to  the  second  mode,  while  the  open  squares  and  dashed  line  corresptmd 
to  the  third.  The  prime  mover  is  filled  with  helium  at  a  mean  pressure  of  170  kPa. 
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3.  Etevelopment  of  a  Non-Boundary  Laver.  Non-Short  Stack  Model  for  the  Power 
Output  of  a  Prime  Mover 

Although  our  porous  medium  approach  to  prime  movers  offers  decent 
agreement  between  measured  and  predicted  quality  factors,  there  is  room  for 
improvement.  Therefore,  we  have  pursued  a  different  technique,  described  at 
length  by  Swift,7. 8  to  our  prime  mover  data.  The  foundation  of  this  technique  is  the 
calculation  of  the  acoustic  power  generated  by  a  differential  element  of  prime  mover 
surface  area.  Although  the  general  technique  is  described  by  Swift,  we  use  an 
expression  for  the  power  output  not  limited  to  either  the  small  boundary  layer  or 
short  stack  approximations.^  This  expression  is  previously  unpublished.  In  the 
calculations  we  take  into  consideration  the  temperature  dependence  of  the 
thermophysical  properties  of  both  the  gas  and  stack  material. 

The  results  of  the  analysis  are  shown  in  Figs.  5-8,  which  correspond  to  the  same 
data  as  presented  in  Figs.  1-4,  respectively.  These  results  represent  some  of  the  first 
comparisons  between  Swift's  theoretical  approach  and  experimental  measurements. 
Comparing  Fig.  1  to  Fig.  4,  it  is  evident  that  the  present  technique  provides  a  much 
better  fit  to  the  170  kPa  data.  However,  the  two  techniques  do  more-or-less  equally 
well  at  376  and  500  kPa  and  for  the  second  and  third  modes  at  170  kPa.  There  is  a 
tendency  for  the  present  technique  to  under  predict  the  attenuation  at  zero 
temperature  difference. 

Our  analysis  provides  good  overall  prediction  of  the  performance  of  a  prime 
mover  below  onset.  However,  the  true  test  is  predicting  the  performance  above 
onset.  There  are  several  significant  differences  between  prime  movers  above  and 
below  onset.  Most  obvious  is  that  the  acoustic  amplitudes  above  onset  are  quite 
large.  Ratios  of  acoustic  pressure  amplitude  to  mean  gas  pressure  of  1  -  10%  are 
common.  Yet,  our  analysis  is  based  on  linear  acoustics.  Beyond  the  theoretical 
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Figure  5.  Graph  of  1  /Q  vs  the  temperature  difference  across  the  prime  mover 
stack.  The  line  represents  the  results  of  the  calculations  based  on  Swift's  theoretical 
approach.  The  data  are  the  same  as  presented  in  Fig.  1;  the  prime  mover  is  filled 
with  helium  at  a  mean  pressure  of  170  kPa. 


Figure  6.  Graph  of  1/Q  vs  the  temperature  difference  across  the  prime  mover 
stack.  The  line  represents  the  results  of  the  calculations  based  on  Swift's  theoretical 
approach.  The  data  are  the  same  as  presented  in  Fig.  2;  the  prime  mover  is  filled 
with  helium  at  a  mean  pressure  of  376  kPa. 
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Figure  7.  Graph  of  1  /Q  vs  the  temperature  difference  across  the  prime  mover 
stack.  The  line  represents  the  results  of  the  calculations  based  on  Swift’s  theoretical 
approach.  The  data  are  the  same  as  presented  in  Fig.  3;  the  prime  mover  is  filled 
with  helium  at  a  mean  pressure  of  500  kPa. 


Figure  8.  Graph  of  1  /Q  vs  the  temperature  difference  across  the  prime  mover 
stack  for  the  second  and  third  longitudinal  modes  showing  the  results  of  the 
calculations  based  on  Swifts  theoretical  approach..  The  x’s  and  solid  line  correspond 
to  the  second  mode,  while  the  open  squares  and  dashed  line  correspond  to  the  third. 
The  data  are  the  same  as  presented  in  Fig.  4;  the  prime  mover  is  filled  with  helium 
at  a  mean  pressure  of  170  kPa. 
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complications  added  by  introducing  nonlinear  effects,  there  are  a  number  of 
experimental  complications.  For  instance,  the  acoustic  displacement  amplitudes 
approach  the  length  of  the  heat  exchangers.  Also,  because  of  the  l£irge  transport  of 
heat,  there  is  no  assurance  that  the  temperature  gradient  along  the  prime  mover 
stack  will  be  uniform. 

The  next  logical  step  in  this  research  is  to  investigate  prime  movers  above 
onset,  but  still  in  the  linear  acoustics  regime.  This  study  is  discussed  in  the  next 
section. 

4.  Study  of  the  Transition  to  Steady  State  of  a  Prime  Mover  Above  Onset 

The  natural  progression  of  our  investigations  of  prime  movers  below  the  onset 
of  self-oscillation  is  to  study  the  transition  to  steady  state  oscillation  above  onset. 
This  work  has  been  initiated  in  FY  1991.  The  below  onset  work  consists  of 
determining  the  net  attenuation  coefficient  of  a  prime  mover  as  a  function  of 
applied  temperature  gradient.  The  attenuation  coefficient  is  determined  from  the 
frequency  response.  The  goal  of  the  above  onset  work  is,  also,  to  determine  the  net 
attenuation  coefficient  as  a  function  of  applied  temperature  gradient.  Above  onset, 
the  attenuation  (gain?)  coefficient  can  be  determined  from  the  rise  time  of  the 
acoustic  oscillations.  However,  above  onset  the  rise  time  is  time  dependent.  It  is 
small  until  some  mechanism  sets  in  to  establish  steady  state,  at  which  point  it  is 
infinite.  Wheatley,  et  al  made  similar  measurements,^ 0  although  they  were  of  a 
preliminary  nature  and  little  quantitative  analysis  was  performed. 

Our  measurements  are  made  in  a  prime  mover  fitted  with  a  butterfly  valve. 
With  the  valve  closed  the  temperature  gradient  across  the  prime  mover  stack  can  be 
raised  above  the  critical  value  required  for  onset  without  oscillations  starting. 
When  the  valve  is  suddenly  opened,  the  prime  mover  goes  into  self-oscillation  and 
the  acoustic  amplitude  begins  to  build  up.  Finally,  the  build-up  ceases  and  steady 
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state  is  reached. 

We  are  interested  in  three  regions  of  the  build-up  of  the  acoustic  oscillations. 
First,  in  the  early  stages  of  the  build-up  the  amplitudes  are  low.  From  our  past 
work,  we  expect  that  measurements  made  in  this  region  should  agree  fairly  well 
with  predictions  and  the  rise  time  to  be  constant.  The  second  region  of  interest  is 
the  part  of  the  build-up  where  the  rise  time  begins  to  increase,  deviating  from  its 
initial  value.  Knowing  at  what  acoustic  amplitude  the  rise  time  starts  to  deviate 
from  its  initial  value  should  give  some  insight  into  what  mechanisms  cause  this 
deviation.  Finally,  we  are  interested  in  the  energy  balance  at  steady  state.  In  other 
words,  we  want  to  make  sure  that  we  can  account  for  the  dissipation  of  all  the 
energy  supplied  by  the  prime  mover.  At  this  point  we  are  investigating  the  low 
amplitude  region.^  ^  The  rise  time  of  the  oscillations  is  measured  as  a  function  of 
the  temperature  gradient  and  mean  gas  pressure  in  the  prime  mover.  The  results 
are  analyzing  using  the  technique  based  on  Swift’s  approach  discussed  in  Section  2. 

We  have  found  that  soon  after  the  oscillations  begin  harmonics  start  to  appear 
in  the  spectrum  of  the  waveform  and  the  rise  time  deviates  from  its  inital  value. 
This  behavior  is  indicated  in  Fig.  9  which  shows  a  plot  of  the  envelope  of  the  signal 
amplitude  (solid  line)  as  a  function  of  time.  In  this  data  record,  harmonics  began  to 
appear  in  the  spectrum  after  approximately  0.4  s.  The  dashed  line  is  a  least  squares 
fit  to  the  first  0.4  s  of  the  data  record.  The  two  lines  deviate  significantly  in  the  later 
portions  of  the  record.  Figure  10  shows  the  logarithm  of  the  amplitude  as  a  fimction 
of  time.  As  a  result  of  the  observations  we  have  limited  our  initial  investigations  to 
measuring  the  initial  rise  time  of  the  oscillations,  before  the  spectrum  shows  any 
significant  harmonic  content.  Aside  from  this,  these  observations  have  taught  us 
that  harmonic  generation  plays  a  significant  role  in  the  transition  to  steady  state. 

Results  of  rise  time  measurements  for  a  helium  filled  prime  mover  are  shown 
in  Figs.  11-14.  To  be  consistent  with  the  presentation  of  the  data  acquired  below 
onset,  the  rise  time  has  been  converted  to  1/Q.  Three  points  of  comparison  are 
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Figure  9.  Graph  of  the  envelope  of  the  signal  amplitude  (solid  line)  as  a  function 
of  time.  The  dashed  line  is  a  least  squares  fit  to  the  first  0.4  s  of  the  data  record. 


Figure  10.  Graph  showing  the  logarithm  of  the  envelope  of  the  signal  amplitude 
as  a  function  of  Hme. 
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Figure  11.  Data  obtained  for  238  kPa  plotted  as  a  function  of  the  temperature 
difference  across  the  stack.  The  solid  line  is  the  theoretical  result. 


Figure  12.  Data  obtained  for  308  kPa  plotted  as  a  function  of  the  temperature 
difference  across  the  stack.  The  solid  line  is  the  theoretical  result. 


13 


p 


Figure  13.  Data  obtained  for  376  kPa  plotted  as  a  function  of  the  temperature 
difference  across  the  stack.  The  solid  line  is  the  theoretical  result. 


H^re  14.  Data  obtained  for  500  kPa  plotted  as  a  function  of  the  temperature 
difference  across  the  stack.  The  solid  line  is  the  theoretical  result. 
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evident.  First,  considering  that  there  are  no  adjustable  parameters  in  the  model,  the 
overall  agreement  is  good,  especially  at  the  higher  pressures.  Second,  the  agreement 
between  the  predicted  and  measured  onset  temperatures,  where  1/Q  =  0,  for  the 
various  mean  pressures  is  inconsistent.  The  third  point  of  comparison  is  that  the 
slopes  of  the  predicted  and  measured  results  are  quite  close.  This  pK)int  is  illustrated 
in  Figs.  15-18.  The  predicted  and  experimental  results  presented  in  Figs.  11-14  were 
extrapolated  to  find  the  temperature  difference  at  onset.  The  results  are  replotted  in 
Figs.  15-18  as  a  function  of  the  number  of  degrees  above  the  extrapolated  onset.  In 
other  words,  the  data  and  theory  have  been  forced  to  fit  at  onset.  As  is  evident,  the 
agreement  is  very  good  at  the  higher  pressures  but  not  quite  as  good  at  the  lower 
pressures.  Recall  that  the  agreement  below  onset  followed  the  same  trend.  The 
reason  for  this  tendency  is  unknown.  However,  one  possible  cause  has  to  do  with 
imperfections  in  the  stack  construction.  At  the  lowest  pressure,  the  plate  spacing  to 
penetration  depth  ratio  is  on  the  order  of  two  or  less,  which  is  just  at  the  critical 
value.  The  heat  transport  occurs  in  a  region  about  a  penetration  depth  away  from 
each  plate.  Therefore,  imperfections  in  the  plate  spacing  could  alter  the  volume  of 
gas  available  for  heat  transport  significantly  at  these  pressures.  However,  these 
imperfections  would  have  less  effect  at  higher  pressures  where  the  plate  spacing  to 
penetration  depth  ratio  is  larger. 

The  accurate  prediction  of  the  slope  is  more  important  than  an  accurate 
prediction  of  the  onset  temperature.  The  inability  to  predict  the  precise  onset 
temperature  may  be  tied  to  errors  in  predicting  the  attenuation  of  the  entire  prime 
mover  at  zero  temperature  difference.  Whereas  the  ability  to  predict  the  slope  well 
indicates  that  the  model  is  able  to  predict  the  change  in  power  output  with  applied 
temperature  difference  well.  In  other  words  our  understanding  of  the  stack 
performance  and  how  that  performance  changes  with  changing  temperature 
gradient  is  good.  We  are  currently  taking  data  with  the  prime  mover  filled  with 
argon  in  attempt  to  further  investigate  the  nature  of  the  discrepancies. 
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Figure  15.  Data  obtained  for  238  kPa  plotted  as  a  funcHon  of  the  number  of  degrees 
above  the  extrapolated  onset. 


Figure  16.  Data  obtained  for  308  kPa  plotted  as  a  function  of  the  number  of  degrees 
alwve  the  extrapdated  onset. 
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Figure  17.  Data  obtained  for  376  kPa  plotted  as  a  function  of  the  number  of  degrees 
above  the  extrapolated  onset. 


Figure  18.  Data  obtained  for  500  kPa  plotted  as  a  function  of  the  number  of  degrees 
alwve  the  extrapolated  onset. 


5.  Study  of  Energy  Distribution  and  Dissiparion  in  Finite  Amplitude  Standing 
Waves 


In  addition  to  the  harmonic  generation  discussed  in  Section  4,  we  have 
previously  observed  nonlinear  effects  in  our  prime  mover  above  onset.'2,i3  These 
effects  are  manifested  in  the  generation  of  harmonics  of  the  fundamental  frequency, 
thus  taking  energy  out  of  the  fundamental  mode.  Because  typical  standing  wave 
thermoacoustic  engines  are  designed  to  operate  optimally  at  the  fundamental 
frequency,  the  transfer  of  energy  to  higher  harmonics  represents  a  parasitic  loss. 
Some  questions  arise.  How  much  of  the  input  energy  is  lost  from  the  fundamental 
mode?  Is  there  a  way  to  prevent  this  loss?  Some  of  these  questions  can  be  answered 
by  understanding  the  steady  state  energy  balance  in  the  prime  mover.  As  a  first  step 
toward  understanding  energy  balance  in  the  prime  mover,  we  have  begun  an 
investigation  of  the  energy  balance  in  a  simpler  system  -  a  rigid  tube. 

Our  efforts  in  this  area  began  with  the  observation  that  the  waveforms 
observed  in  our  prime  mover  are  very  similar  to  those  observed  by  Coppens  and 
Sanders.^'*  Their  system  consisted  of  a  rigid  tube  driven  with  a  piston  source  at  one 
end.  We  have  created  a  similar  setup  (using  their  origf.nal  tube  and  source)  and 
made  measurements  of  the  steady  state  energy  distribution  in  finite  amplitude 
standing  waves. 

The  first  measurements  were  made  on  an  empty  tube.  The  tube  was  driven  at 
its  fundamental  resonance  frequency.  The  steady  state  power  input  to  the  tube  by 
the  piston  was  computed  from  measurements  of  the  acceleration  of  the  piston  and 
the  acoustic  pressure  amplitude  at  the  piston  location.  This  power  was  compared  to 
the  steady  state  acoustic  power  dbsipated  by  thermal  and  viscous  losses  at  the  tube 
walls  by  the  fundamental  and  higher  harmonics.  The  results  are  shown  in  Figs.  19 
and  20.  In  these  figures,  ^  represents  the  power  dissipated  by  the  nf^  harmonic 
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(where  n  =  1  corresponds  to  the  fundamental).  Einput  is  the  total  power  input  to  the 
tube  by  the  piston.  It  can  be  seen  that  at  the  highest  drive  levels,  only  85%  of  the 
power  remains  in  the  fundamental.  Conversely,  15%  of  the  power  delivered  at  the 
fundamental  is  dissipated  by  higher  harmonics,  the  result  of  nonlinear  effects. 
Figure  21  shows  the  ratio  of  the  total  dissipated  power  E=Ei  +E2  +  -  -  to  Empm.  The 
fact  that  this  ratio  is  very  close  to  unity  implies  that  our  system  behaves  as  expected. 
Figure  22  shows  a  typical  waveform  and  spectrum  of  the  acoustic  pressure  at  the 
highest  drive  levels.  Notice  the  severe  distortion  of  the  waveform  and  the 
associated  rich  spectral  content. 

The  next  set  of  measurements  was  made  on  an  obstructed  tube.  The 
obstruction,  shown  in  Fig.  23,  is  an  insert  located  at  the  midpoint  of  the  original 
empty  tube.  As  seen  from  Table  1,  the  obstruction  changes  the  resonance 
frequencies  of  the  tube  so  that  the  harmonics  of  the  fundamental  (drive)  frequency 
no  longer  match  the  overtones  of  the  tube.  It  detunes  the  resonator.  The  same 
measurements  were  repeated  as  with  the  empty  tube.  Figure  24  shows  a  typical 
waveform  and  spectrum  at  the  same  drive  levels  as  shown  in  Fig.  22.  (Note  that  the 
time  and  frequency  spans  are  different  from  those  in  Fig.  22.)  The  waveform  is 
much  less  distorted.  This  can  be  seen  visually  or  by  the  relatively  few  numbers  of 
harmonics  in  the  spectrum.  Clearly,  the  detuning  cleans  up  the  waveform.  But 
what  effect  does  this  have  on  the  energy  distribution?  Do  we  still  know  where  all 
the  energy  is  going? 

The  second  harmonic  of  the  pressure  waveform  is  more  than  35  dB  below  the 
fundamental  amplitude  and  so  accounts  for  less  than  a  few  hundredths  of  a  percent 
of  the  total  acoustic  energy.  Therefore,  E  can  be  approximated  by  Ei.  Figure  24  shows 
E/Einpat  for  the  obstructed  tube.  There  is  a  systematic  decrease  in  the  ratio  from 
approximately  unity  to  0.93  as  the  drive  amplitude  increases.  For  the  reasons  given 
at  the  beginning  of  the  paragraph,  this  energy  does  not  show  up  in  other  acoustic 
modes.  Where  the  unaccounted  for  energy  goes  is  not  known.  There  may  be 
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Figure  19.  Measured  ratio  of  Ei/Einpqi  as  a  function  of  the  piston  Mach  nuntber  for 
the  empty  tube. 
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Figure  20.  Measured  ratio  of  E2+E3+  .7Eiiii«ii  as  a  function  of  the  piston  Mach 
number  for  the  empty  tube. 
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Figure  22,  A  typical  waveform  (top)  and  spectrum  (bottom)  of  the  acoustic  pressure 
at  the  highest  drive  levels  in  the  empty  tul^. 
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Figure  23.  Illustration  of  the  obstruction  that  is  inserted  into  the  empty  tube. 


Mode 

Empty  Tube 

Obstructed  Tube 

fn  (Hz)  fn/fi 

fn  (Hz)  fn/fi 

1 

238 

1.000 

227 

1.000 

2 

477 

2.004 

487 

2.144 

3 

715 

3.004 

686 

3.019 

4 

955 

4.013 

974 

4.289 

5 

1194 

5.017 

1153 

5.077 

Table  1.  Measured  resonance  frequencies  for  the  first  five  modes  of  the  empty  and 
obstructed  tubes. 
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turbulence  generated  at  the  tube /obstruction  transitions.  Also,  there  are  cavities 
between  the  inner  tube  wall  and  the  outer  obstruction  wall  that  may  be  contributing 
to  the  problem.  We  are  in  the  process  of  redesigning  the  resonator  to  eliminate 
these  cavities. 

In  spite  of  the  systematic  decrease  in  the  ratio,  the  results  with  the  obstructed 
tube  clearly  are  better  than  the  15%  energy  loss  from  the  empty  tube.  The  results 
suggest  that  resonators  with  nearly  harmonic  overtones  are  to  be  avoided  in 
thermoacoustics  work.  This  conclusion  may  seem  painfully  obvious,  but  now  we 
have  some  ideas  of  the  magnitude  of  the  problem.  There  is  also  probably  something 
to  learn  from  the  systematic  decrease  in  E for  the  obstructed  tube.  Finally,  we 
are  able  to  account  for  the  energy  balance  in  finite  amplitude  standing  waves  (in 
empty  tubes  at  least).  This  ability  will  be  applied  to  steady  state  oscillations  in  prime 
movers  in  FY  1992. 
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Figure  25.  Measured  ratio  of  E/Einpot  as  a  function  of  the  piston  Mach  number  for 
the  obstructed  tube. 
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